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AUTHOR'S  SUMMARY 


The  separated  flow  of  an  ideal  fluid  around  delta  wings  of  low  aspect  ratio 
with  rounded  leading  edges  is  calculated  by  an  approximation  to  the  theory  of 
slender  bodies.  The  numerical  method  employed  is  based  upon  a  mathematical  model 
of  inviscid  separation  from  a  smooth  surface.  Results  of  calculations  are  given 
for  symmetrical  flow  around  wings  whose  cross-sections  are  ellipses.  Particular 
attention  is  paid  to  the  investigation  of  the  influence  of  the  position  of  the 
line  of  separation  in  the  vicinity  of  the  rounded  edge  of  the  wing  on  the  con¬ 


figuration  of 
the  wing. 


the  vortex  sheet  and  to  the  overall  aerodynamic  characteristics  of 
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Experimental  results  show  that  the  characteristics  of  separated  flow  around 

delta  wings  of  low  aspect  ratio  depend  substantially  upon  the  form  of  the  leading 

1  .23 

edges  .  Existing  methods  of  calculation  of  such  flows  ’  within  the  constraints 

of  an  ideal  fluid  are  based  upon  the  assumption  of  sharp  edges  on  the  wing,  even 

when  the  thickness  of  the  wing  is  taken  into  consideration^.  Although  such  a 

constraint  on  the  shape  of  the  wings  does  allow  a  substantial  simplification  of 

the  task  (complying  with  the  Joukowski  hypothesis  for  sharp  wings  ensures  a 

unique  solution),  the  area  of  application  of  these  methods  is  restricted.  In 

instances  when  these  methods  are  applied  to  the  calculation  of  flow  around  delta 

wings  with  rounded  leading  edges,  they  do  not  achieve  solutions  of  the  required 

accuracy.  Rounding  of  the  wings  gives  rise  to  a  finite  range  of  angles  of  attack 

in  which  the  flow  is  generally  unseparated.  At  large  angles  of  attack,  separated 

flow  occurs,  but  the  supplementary  increment  in  lifting  force  is  significantly 

less  than  in  the  case  of  flow  around  the  same  wings  with  sharp  edges*. 

For  the  calculation  of  the  characteristics  of  separated  flow  around  wings 
of  low  aspect  ratio,  for  which  the  radii  of  rounding  of  the  edges  is  much  greater 
than  the  thickness  of  the  boundary  layers,  a  method  may  be  used,  based  upon  a 
moael  of  inviscid  separated  flow  around  bodies  with  a  smooth  surface.  This  model 
is  characterised  by  the  non-unique  solutions  obtained.  In  particular,  there  may 
occur  a  single-parameter  set  of  solutions  with  the  position  of  the  lines  of 
separation  as  an  independent  parameter,  as,  for  example,  in  the  problem  of 
symmetrical  separated  flow  around  a  thin  circular  cone^. 

The  aim  of  this  article  is  the  construction  of  single-parameter  sets  of 
solutions  to  the  problem  of  symmetrical  separated  flow  around  delta  wings  of  low 
aspect  ratio,  whose  cross-sections  are  ellipses  of  small  relative  thickness;  also 
an  analysis  on  the  basis  of  these  solutions  of  the  influence  of  the  position  of 
the  line  of  separation,  the  rounding  of  the  edges  and  the  thickness  of  the  wing  on 
the  vortex  system.  The  aim  is  also  to  obtain  the  overall  aerodynamic 
characteristics. 

The  question  of  the  extraction  of  the  unique  solution  from  the  correspond¬ 
ing  single-parameter  set  of  solutions  goes  beyond  the  limits  of  the  present 
investigation.  The  answer  to  this  may  be  found  either  experimentally  (the  deter¬ 
mination  of  the  position  of  the  line  of  separation  in  each  particular  case),  or 
theoretically  by  a  special  investigation  of  both  the  regions  of  flow  in  which 

viscous  forces  play  a  substantial  part,  and  the  application  of  the  corresponding 

.  .  .  .  .  6 

procedures  for  bringing  together  the  solutions  in  these  regions  . 
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1  Sfan'np  the  problem.  Assume  a  right-angled  Cartesian  system  of  co-ordinates 

with  its  origin  at  the  wing  apex.  Let  the  axis  Oz  be  directed  along  the  axis  of 
symmetry  of  the  wing,  the  axis  Ox  along  the  span  and  the  axis  Oy  perpendicular 
to  the  axes  Oz  and  Ox  .  The  wing  is  subject  to  a  flow  of  gas  around  it  and  at 
an  angle  of  attack  a  with  a  velocity  vector  at  infinity,  and  lying  in  the 

(y,z)  plane  (Fig  1). 


In  accordance  with  the  results  of  the  Ward-Jones  theory  of  slender  bodies  , 
the  velocity  potential  (f)(x,y,z)  at  distances  of  the  order  of  the  typical  dimen¬ 
sion  of  the  cross-section  may  be  expressed  in  the  following  form: 


<t)(x,y,z) 


+  <j)j  (x,y,z) 


where  with  M  <  1 


2TTg  (z)  =  s'(z)  In - 2 


and  with  M  >  I  ,  S'(0)  =  0 


-  1 

2TTgQ(z)  =  S'(z)  In  - 2 - 


1 

0 

j  S"(u)  sign(u  -  z)  In  |p 
0 


z 

j  S"(y)  In  (z  -  u)du  . 
0 


z  I  dp 
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Here  is  the  quantity  M  of  the  airstream,  S(z)  the  area  of  the  wii 

cross-section,  0  ^  z  <  1  .  The  function  (p  is  harmonic  in  the  variables 

2  2 

x,y($V  4>V  =  0)  and  for  x  +  y  may  be  expressed  in  the  form: 

Ixx  lyy  ^ 

'^1  ^  .  S'(z)  ,  ,r2  2 

—  =5:  ay  +  — -  InVx  +  y 


/ 

Also,  the  expression  11 - —Jz  +  (()j(;<,y,z)  satisfies  the  condition  for 

absence  of  flow  through  the  surface  of  the  wing  and  the  boundary  conditions  on  the 
vortex  sheets. 


2  +  (()j(;<,y,z)  satisfies  the  condition  for 


After  introducing  the  time-variable  t  =  z/w  ,  the  three-dimensional 

oo 

stationary  problem  of  the  separated  flow  around  the  wing  under  examination  reverts 
to  the  two-dimensional  unsteady  (similarity)  problem  of  the  separated  flow  of  an 
incompressible  fluid  which  expands  uniformly  in  time  around  the  elliptical  contour 
Lq  (Fig  2a). 


Fig  2 

For  a  solution  of  the  above  we  use  the  methods  of  the  theory  of  functions 
of  a  complex  variable.  In  the  plane  of  the  cross-section  of  the  wing  we  introduce 
the  complex  variable  C  =  x  +  iy  and  examine  the  flow  picture  at  an  arbitrary 
fixed  moment  in  time.  Let  the  configuration  and  intensity  of  the  vortex  sheets 
Lj  and  be  given  in  the  parametric  form 


4],2(F,t) 


0  <  Irl 


where  F  is  the  circulation  of  the  section  of  the  vortex  sheet,  measured  from  the 
core  (the  actual  parameter).  Their  magnitudes  or  signs  at  the  points  of  separa¬ 
tion  of  the  flow  are  denoted  both  here  and  in  the  following  text  by  an  asterisk  and 
subscript . 
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Let  us  formulate  the  folloviing  problem:  it  is  required  to  find  the  complex 
velocity  of  flow  V(c»t)  which  is  regular  outside  the  contour  L  =  *  ^2  ' 

and  possessing  the  specified  discontinuity  on  the  vortex  sheets  and 

[V^(?,t)  -  V_(C,t)i  =  Y,  " 

where  "+"  and  indicate  the  signs  of  V  from  the  two  sides  of  the  sheet, 

Y  represents  the  linear  density  of  the  circulation;  the  complex  velocity  is  also 
to  satisfy  the  boundary  conditions  for  the  absence  of  flow  on  the  contour 
and  the  attenuation  of  perturbations  at  infinity 


V(^,t) 


-  law 


The  sought-after  complex  velocity  satisfying  all  the  stated  conditions  may 
be  readily  obtained  if  the  corresponding  complex  potential  of  the  flow  is  known. 
To  find  the  latter  we  use  the  conformal  mapping  of  the  outside  of  the  ellipse 
with  semi-major  axes  a(t)  and  b(t)  in  the  plane  of  the  variable  ;  onto  the 
outside  of  the  circle  of  radius  R(t)  in  the  plane  of  the  complex  variable 
0  =  C  +  in  (Fig  2b) 

0  =  iU  *  -  d^)-]  . 


where  d  -  a/l  -  ,  R  =  a(l  +  e)/2  ,  t  =  b/a  ,  the  relative  thickness  of  the 

wing.  The  inverse  mapping  has  the  form 


=  a  + 


4o 


In  the  particular  case  of  a  symmetrical  scheme  of  separated  flow  we  have  for 
the  complex  potential  W(a,t)  of  the  flow  under  investigation 


da 


tt  1  •  I  E\ 

W  *  ca  —  In  a  -  law  |  o - 

dt 


(a  -  Oq) (a  +  - 


(a  +  Oq)  ^ 
'^0 


2s 


Both  here  and  in  the  following  text,  the  variable  with  the  bar  indicates 
the  complex  conjugate  of  the  quantity.  The  link  between  the  complex  potential 
and  the  required  flow  velocity  in  the  plane  of  the  variable  c  is  expressed  in 
the  usual  manner 
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I 


V(c,t) 


(1.2) 


As  soon  as  the  complex  flow  velocity  for  the  arbitrary  fixed  moment  in  time 

is  obtained,  the  configuration  and  intensity  of  the  vortex  sheets  at  subsequent 

moments  of  time  may  be  obtained  as  a  result  of  Cauchy's  solution  for  singular 

8 

integro-dif ferentiai  equations 


^^1  2 

at  vl^i ,t]  , 


0  <  |ri  < 


which  are  able  to  fulfil  the  conditions  of  continuity  of  pressure  and  normal 
component  of  flow  velocity  on  the  vortex  sheets. 

2  Method  of  solution.  The  Cauchy  solution  as  expressed  above  ( 1 .) )  to  (1.3)  is 

applied  to  the  vortex  sheets  at  a  given  fixed  moment  of  time  in  the  flow.  To 

obtain  the  sections  of  the  vortex  sheets  leaving  the  surface  of  the  wing  and 

entering  the  stream  in  the  interval  of  time  from  the  given  value  to  the  next,  it 
9 

was  proposed  to  use  the  equation  for  the  density  of  the  circulation  at  the  point 
of  separation  y  and  the  expression  for  the  velocity  of  shift  of  the  point  of 
separation  v*  along  the  smooth  contour  of  a  body  under  the  condition  p  =  r  , 
where  o  is  the  radius  of  curvature  of  the  vortex  sheet,  and  r  is  the  radius 
of  curvature  of  the  surface  of  the  body. 


(2.1) 


Here  v  ,  v  are  the  tangential  components  of  velocity  from  the  outer 
s-*-  s-  " 

and  inner  sides  of  the  sheet  respectively,  v  is  the  normal  component  of  velocity 

n 

on  the  sheet,  and  s  is  the  length  of  the  arc. 

As  shown^’'*^,  the  integration  of  the  equation  (2.1)  with  p  =  r  in  common 
with  equations  (1.1)  to  (1.3)  leads  to  the  achievement  of  a  unique  solution  satisfy¬ 
ing  the  Brillouin-Villat  conditions:  at  the  points  of  separation  there  is  zero 
gradient  of  pressure  and  the  curvature  of  the  sheet  equals  the  curvature  of  the 
body.  Analysis  of  the  asymptote  of  the  flow  in  the  vicinity  of  the  point  of 
separation  reveals  the  indeterminacy  of  the  expression  in  curved  brackets  in  equa¬ 
tion  (2.1),  arising  in  the  case  when  solutions  are  sought  which  do  not  satisfy 
the  Brillouin-Villat  conditions.  As  is  known^,  these  solutions  are  characterised 
by  an  infinite  curvature  of  the  vortex  sheet  at  the  point  of  separation,  and  con¬ 
sequently  by  an  infinite  pressure  gradient  at  this  point. 
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The  discovery  of  the  indeterminacy  in  (2.1)  leads  to  an  equation  in  which  an 
arbitrary  function  occurs,  and  in  the  case  of  self-modelling  flow,  an  arbitrary 
constant*^,  indeterminable  within  the  constraints  of  an  ideal  fluid.  In  deter¬ 
mining  the  significance  of  this  constant,  solutions  may  be  obtained  for  similarity 
problems  of  separated  flow  with  a  different  position  of  the  separation  points. 

Such  an  implicit  operation  for  obtaining  the  positions  of  separation  was  used^  to 
obtain  the  single-parameter  set  of  solutions  for  the  inviscid  symmetrical  separ¬ 
ated  flow  around  a  thin  circular  cone. 

In  the  present  work,  an  explicit  programme  for  the  separation  points  on  an 
elliptical  contour  was  used.  Even  if  such  an  approach  leads  to  a  disturbance  in 
the  equations  of  motion  in  the  process  of  determining  the  solution,  it  is  justi¬ 
fied  for  self-modelling  problems,  characterised  by  the  conditions 

=  0  ,  V*  =  0  ,  V*  =  ('"'g^.) 

The  process  of  determining  the  solutions  in  this  case  is  accelerated,  inso¬ 
far  as  there  is  an  absence  of  the  oscillatory  movement  of  the  points  of  separation 
along  the  contour  of  the  body,  which  takes  place  while  obtaining  a  solution  with 
an  implicit  programme  of  the  position  of  separation. 

The  impossibility  of  obtaining  a  solution  precisely  satisfying  the 
Bri llouin-Villat  conditions  is  in  some  way  a  shortcoming  of  the  approach  under 
investigation.  However,  the  position  of  the  point  of  separation  corresponding  to 
the  specified  condition  may  be  determined  with  sufficient  accuracy  on  account  of 
the  absence  of  convergence  in  the  process  of  determining  the  earlier  (upstream) 
positions  of  the  separation  points. 

In  other  respects  the  numerical  method  used  corresponds  to  the  one  quoted^, 
including:  the  discretization  of  the  vortex  sheet;  the  integration  of  the  system 
of  equations  of  motion  of  the  discrete  vortices  according  to  a  regularised  scheme 
having  a  first  order  of  accuracy;  the  use  of  a  model  of  the  core  of  the  sheet 
(vortex  and  cut);  the  improved  approximation  of  velocity  along  the  normal  to  the 
vortex  sheet,  at  intervals  less  than  the  steps  of  the  discretization. 

3  Results  of  the  calculations.  In  view  of  the  large  number  of  parameters  in 
the  problem  under  investigation,  calculations  were  provided  at  fixed  values  of 
the  relative  angle  of  attack  a  *  o/x  =  1-5  ,  where  x  is  the  half-angle  at  the 
apex  of  the  wing  in  the  plane  (x,z)  for  three  wings  with  relative  thicknesses 
e  =  1/7;  I/IO;  1/14.  In  accordance  with  the  experimental  data*',  with  a  fixed 
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value  of  che  relative  angle  of  attack  the  lines  of  separation  on  elliptical  cones 
with  a  diirinished  relative  thickness  are  directed  towards  a  generator  passing 
through  the  tips  of  the  ellipses,  and  ending  on  the  upper  side  of  their  surface. 

Arising  from  this,  the  vicinity  of  the  ends  of  the  large  semi-major  axes 
of  the  ellipses  was  chosen  to  represent  the  range  of  variation  in  the  position 
of  the  line  of  separation.  Let  us  agree  in  future  to  refer  to  this  as  the  vicin¬ 
ity  of  the  rounded  edge,  and  the  geometrical  position  of  the  tips  of  the  ellipses, 
at  which  the  curvature  is  a  maximum,  as  the  edge  itself. 

Let  us  introduce  into  the  section  z  =  const  a  curvilinear  system  of 
co-ordinates  x^,  y|  .  The  abscissa  x^  is  referred  to  the  edge  along  the  con¬ 
tour  of  the  ellipse,  and  the  ordinate  y^  to  the  outer  normal  to  the  contour. 

We  will  designate  the  co-ordinates  of  the  point  of  separation  in  this  system  of 
co-ordinates  by  x*,  y|  —  0  .  For  the  configuration  of  the  vortex  sheet  in  the 
vicinity  of  the  point  of  separation,  the  approximation  is  justifiably  made 


X.) 

I 


3/2 


X]  <  x*  <  X** 


and  correspondingly  in  the  limiting  case  of  separated  flows  satisfying  the 
Brillouin-Villat  conditions. 


Vj  5:  Cj(x**  -  X|)^''“,  Xj  <  X**  , 

where  x**  ,  y**  =  0  are  the  co-ordinates  of  the  point  of  separation  in  this 

limiting  case. 

The  coefficients  ‘^q  »  j  represent  functions  of  a  ,  e  ,  and  the 
coefficient  also  represents  a  function  of  (x**  -  x*)  . 

At  the  right  of  Fig  3  the  vicinities  of  the  rounded  edges  are  shown  together 
with  those  regions  under  investigation  of  the  positions  of  the  lines  of  separation 
on  all  three  wings.  The  upper  boundaries  of  the  regions  (see  Fig  3)  are  arbi¬ 
trarily  chosen. 
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with  regard  Co  the  lower  boundaries  of  the  region  under  investigation  it  is  pos¬ 
sible  almost  to  take  them  beyond  the  limits  of  existence  of  solutions  (they 
correspond  closely  to  solutions  fulfilling  the  Brillouin-Villat  conditions  at  the 
point  of  separation).  As  already  noted  above,  the  accurate  achievement  of  solu¬ 
tions  satisfying  these  conditions  is  not  possible  by  the  method  employed  in  this 
work  of  directly  programming  the  positions  of  the  separation  points.  The 
behaviour  of  the  lower  boundaries  of  the  regions  of  the  position  of  the  line  of 
separation  under  investigation  is  consistent  with  the  tendency  of  x**  towards 

the  edge  of  the  wing  with  a  decrease  in  relative  thickness  e  and  consequently 

2 

in  the  radius  of  the  rounded  edge  ("e  )  . 

At  the  left  of  Fig  3,  the  change  of  position  of  the  cores  of  the  vortex 
sheets  for  the  ranges  of  positions  of  the  line  of  separation  under  investigation 
is  shown  in  x,y  co-ordinates  (see  Fig  2a)  in  continuous  lines  to  the  same  scale 
The  arrow  indicates  the  direction  of  the  change  in  position  of  the  cores  of  the 
sheets,  this  direction  corresponding  to  the  change  in  the  line  of  separation  from 
the  highest  to  the  lowest  boundaries  of  the  ranges  under  investigation.  The 
symbols  on  the  curves  indicate  the  position  of  the  cores  of  the  vortex  sheets  at 
separation  from  the  edges  (x*  =  0''  .  The  small  displacement  of  the  separation 
points  in  the  vicinity  of  the  rounded  edges  (especially  on  the  upper  surface  of 
the  wing)  leads  to  a  subscan..ial  displacement  of  the  cores  of  the  vortex  sheets, 
the  larger,  the  smaller  the  value  of  e  .  For  comparison,  in  the  same  place  the 

position  is  shown  of  the  core  of  the  vortex  sheet  for  a  triangular  lamina  with 

3  ... 

sharp  edges  for  the  case  when  the  Joukowski  hypothesis  is  fulfilled  at  these 

edges.  The  form  of  distribution  of  the  curves  for  the  various  values  of  c  , 

4 

also  the  analysis  of  the  results  concerning  the  influence  of  small  thickness  on 
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the  characteristics  of  separated  flow  around  delta  wings  with  sharp  leading  edges 
allow  one  to  propose  the  existence  of  a  simple  mechanism  for  the  displacement  of 
the  vortex  spiral  along  the  y  ordinate  on  the  upper  surface  of  the  wing.  If 
ve  subtract  from  the  ordinate  of  the  position  of  the  core  of  the  vortex  sheet  the 
ordinate  of  the  upper  surface  of  the  wing  with  the  same  abscissa,  then  we  obtain 
the  curves  shown  at  the  lower  left  of  Fig  3.  The  trend  of  the  curves  towards  the 
'plane'  (y  =  0)  of  the  upper  surface  of  the  wing  is  consistent  with  the  tendency 
of  those  solutions  obtained  for  the  position  of  the  cores  of  the  vortex  sheets  to 
approach  that  for  a  plane  triangular  lamina  with  sharp  leading  edges,  with 
reduced  relative  thickness  and  under  the  pondition  of  maximum  early  separation 
from  the  windward  surface  of  the  wing  (meeting  the  Brillouin-Villat  conditions 
at  the  separation  points) . 


In  Fig  4  is  shown  the  dependency  of  the  coefficient  of  normal  force  C  of 

N 

the  wings  under  examination  on  the  co-ordinate  of  the  separation  point  in  the 
vicinity  of  the  rounded  edge.  The  coefficient  of  normal  force  is  related  to  the 
coefficient  of  normal  force  of  a  triangular  lamina  with  sharp  edges: 

CNq/x  =  18.82  with  a  =  1 .5  (see  Ref  3).  The  maximum  value  of  the  coefficient 
of  normal  force,  as  in  the  case  of  separated  flow  around  a  circular  cone^,  occurs 
when  the  Brillouin-Villat  conditions  are  met  at  the  points  of  separation.  With  a 
reduction  in  the  relative  thickness  of  the  wing  e  ,  an  increment  is  observed  in 
this  maximum  value  and  its  tendency  towards  the  coefficient  of  normal  force  for  a 
triangular  lamina  with  sharp  edges. 

Fig  5  illustrates  the  configurations  of  the  vortex  sheets  for  the  extreme 
points  of  the  investigated  ranges  of  positions  of  the  line  of  separation  on  a  wing 
having  relative  thickness  e  =  1/10  . 
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The  investigations  conducted  show  that  even  an  insignificant  displacement 
of  the  position  of  the  lines  of  separation  in  the  vicinity  of  rounded  leading 
edges  on  the  upper  surface  of  the  wing  (in  calculations,  the  displacement  did  not 
exceed  3%  of  the  half-span)  leads  to  a  substantial  reduction  in  the  supplementary 
normal  force  brought  about  by  the  separated  flow.  For  example,  for  a  wing  with 
relative  thickness  e  =  1/10  ,  the  supplementary  normal  force  is  decreased  by  62%. 


The  above  discussion  supports  the  possibility  of  setting-up  a  numerical 
method  for  the  calculation  of  the  characteristics  of  separated  flow  around  wings 
of  small  aspect  ratio,  having  thickness  and  rounded  leading  edges,  within  the 
constraints  of  a  model  of  an  ideal  fluid.  For  its  effective  use,  it  is  necessary 
to  have  means  of  determining  with  the  highest  possible  accuracy  (experimentally 
or  theoretically)  the  position  of  the  lines  of  separation  in  the  vicinity  of  the 
leading  edges  under  given  specific  conditions  of  flow. 
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